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Abstract 

We consider flavor breaking in the the octet and decuplet baryon masses at 
second order in large- chiral perturbation theory, where is the number of 
QCD colors. We assume that ~ ^/Np ~ m^/A ^> m^i^^/A, oeM) where 
Np is the number of light quark flavors, and m^^fi^g/K are the parameters 
controlling SU{Np) flavor breaking in chiral perturbation theory. We con- 
sistently include non-analytic contributions to the baryon masses at orders 
rriq'^ ^ nig Inm^, and {niq lnmq)/N. The m^^^ corrections are small for the re- 
lations that follow from SU{Np) symmetry alone, but the corrections to the 
large- A^ relations are large and have the wrong sign. Chiral power-counting 
and large- A^ consistency allow a 2-loop contribution at order rn^lnruq, and 
a non-trivial explicit calculation is required to show that this contribution 
vanishes. At second order in the expansion, there are eight relations that are 
non-trivial consequences of the 1/A^ expansion, all of which are well satisfied 
within the experimental errors. The average deviation at this order is 7 MeV 
for the AI = mass differences and 0.35 MeV for the AI ^ mass differences, 
consistent with the expectation that the error is of order 1/A^^ ~ 10%. 
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I. INTRODUCTION 



In this paper, we analyze the octet and decuplet baryon masses to second order in a 
simultaneous chiral and expansion, where N is the number of QCD colors. The 
expansion for baryons has a rich structure and significant predictive power: at leading order 
in the expansion, static baryon matrix elements satisfy SU{2Nf) spin-flavor relations 
(where Np is the number of flavors) and the corrections to these relations are highly 
constrained In order to compare the predictions of the 1/N expansion with experiment 

at subleading order, we must consider both corrections and explicit breaking of SU{Np) 
flavor symmetry due to quark masses and electromagnetism. We will use the formalism of 
refs. [^J^ that takes Np ^ N ^ 1 and makes use of an explicit effective lagrangian that 
keeps the SU{Np) flavor symmetry manifest order- by-order in the expansion. In order 
to determine which terms in the double expansion to keep, we will expand assuming that 

1 1 rrig niu d 

^-^-X»^»«EM, (1.1) 

where A is the chiral expansion parameter and oem is the electromagnetic coupling.^ The 
baryon magnetic moments were analyzed in the same expansion in ref. and found to 
be in excellent agreement with experiment. Different expansions for various baryon 
observables have also been considered in refs. PJ9|-[TI|. In the concluding section we will 



briefly compare our results with those of ref. ||T^ , which also considers baryon masses 



One important feature of the present work is that we consistently include the chiral 
loop contributions in our expansion. The leading non-analytic contributions to the baryon 
mass differences are order m^/^, {mq\mmq)/N ^ and m^lnmg. Power-counting and large- 

consistency arguments allow a 2-loop contribution at order mqlnmq, and a non-trivial 
explicit calculation is required to see that such a contribution does not appear. The order 
^3/2 corrections are calculable, and the result is that they are small for the relations that 
follow from SU{Np) symmetry alone, but the corrections to the large- relations are large 
and have the wrong sign. While this may indicate that the expansion parameters are not 
sufficiently small in nature to believe this expansion, we note that there are higher-order 
effects that are expected to substantially reduce these corrections. Also, the m^/^ corrections 
can be cancelled at higher orders, giving results that agree well with experiment. The 
coefficients of the {rriq In ■mq)/N and rrig In rriq corrections are not calculable in terms of 
measured couplings, and we have included them as arbitrary parameters. 

Our final results include corrections up to order and nis/N for the A J = mass 
differences; and order (m^ — mu)ms, {rrid — mu)/N, and Oem for the A/ ^ mass differ- 
ences. At this order there are eight relations which are non-trivial consequences of the 1/A^ 
expansion. (That is, they do not follow from flavor symmetry considerations alone.) These 
relations agree well with experiment, and the remaining deviations are consistent with the 
expectation that the leading corrections are order l/N"^. 



See section IIC for a discussion of the last inequality. 
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This paper is organized as follows. In section 2, we review the expansion used in this 
paper. In section 3, we present our results, and section 4 contains our conclusions. The 
details of a 2-loop computation are contained in an appendix. 

II. THE EXPANSION 

The expansion for baryons has a good deal of predictive power even at subleading 
orders, but some aspects of the expansion are rather subtle. In this section, we review the 
main ingredients of this expansion. 

A. Baryon quantum numbers 

The expansion makes sense only for baryons with spin J ~ 1, since baryons with 
spin J N have width of order A^. One conceptual subtlety in the 1/N expansion for 
baryons is that even for fixed Np > 3 the number of baryon states with given spin J ~ 1 
grows with N. When Np = 2, the baryons have quantum numbers I = J = and it 

is clear that one should identify the states of lowest spin and isospin with the corresponding 
baryon states at = 3. When Np > 3, the situation is more complicated. The Young 
tableaux for the SU{Np) representation of spin J baryons is shown in fig. 1; for > 3, 
this representation contains many states that may be identified with a given baryon state 
in the world with A^ = 3. Choosing any subset of the baryon states that exist for A^ > 3 to 
represent the baryons at A^ = 3 breaks the flavor symmetry explicitly. The "extra" states 
that appear for A^ > 3 are important for computing chiral loops, since they can appear as 
intermediate states; the contributions of all baryon states are required in order to maintain 
flavor symmetry and large- A^ consistency. 

A simple way to handle this situation was pointed out in refs. [^0. The 1/A^ expansion 
can be carried out without selecting any subset of the baryon states for A^ > 3 by writing the 
matrix elements in terms of few-body operators in a spin-flavor Fock space that describes the 
baryon quantum numbers. In this formalism, the coefficient of an r-body operator is at most 
1/N^~^, and there is a simple classification of the operators that allows us to determine the 
A^- and A^i^-dependence of the matrix elements for arbitrary baryon states [H . For operators 
with the flavor numbers of a product of SU (Np) adjoints, any operator can be written as a 
product of 1-body operators of the form 

{XFj^aLX^rV^ (2.1) 

where a'^ (a) is a creation (annihilation) operator in the spin-flavor Fock space, X is a flavor 
matrix and F is a spin matrix (either 1 or ) , with flavor and Lorentz indices contracted in 
all possible ways. We keep a given operator if its matrix element in any baryon state with 
J ~ 1 (for Np > 3) is as large as the order to which we are working. It is not hard to see 
that the largest matrix elements are 

m<{l if^ = l,r = a^ ^2.2) 
L i\ otherwise. 
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Since the operators have definite SU{Nf) quantum numbers, this procedure keeps the fiavor 
symmetry manifest for arbitrary A^. 



B. Large Np 

In ref. it was shown that the large- counting rules for baryons are unaffected when 
the number of flavors is taken to be large, i.e. Np ~ A^. From the point of view of the quark 
and gluon degrees of freedom, the reason for this is that the large-A^ counting rules arise 
from the suppression of many-quark interactions, and this suppression is unaffected by the 
presence of quark loops. This is to be contrasted with the large-iV predictions for mesons, 
which generally rely on the suppression of quark loops and therefore are valid only in the 
limit Np -C A^. 

We will assume that Np ~ ^ 1. This seems sensible, since N = Np = 3 in nature. 
In order to carry out this expansion, we must decide how to extrapolate flavor breaking 
(the quark masses and charges) to a world with Np > 3. Clearly, there are inflnitely many 
choices, and our approach will be similar in spirit to the way we handle the additional 
baryon states that occur: we set up the expansion to be independent of the details of the 
extrapolation. 

We therefore consider extrapolations with arbitrary numbers of individual quark flavors 
Nu,d,s with Nu + Nd + Ns = Np. There are extrapolations where A^^ ~ Np (g = u, d, s), but 
there are also extrapolations where e.g. A^^ ~ 1 and A^^^^ ~ Np. We will evaluate the matrix 
elements for arbitrary Ng, and keep any operator that is as large as the order to which we 
are working on any baryon state for any extrapolation. Thus for example, we keep terms of 
order l/Ng as well as Nq/N. 



C. The double expansion 

The predictions of the large-A^ limit for baryons can be summarized by stating that 
baryon matrix elements obey SU {2Np) spin-flavor relations in this limit. In order to consider 
corrections to this limit, we must take into account both the fact that A^ is flnite, and the 
fact that the SU{Np) flavor symmetry is broken by quark masses and electromagnetism. 
We therefore carry out a simultaneous expansion in 1/A^ and flavor breaking. In order to do 
this, we must decide where to truncate the expansion in the small parameters 1/A^ ~ 1/Np, 
mu,d,s/^, and q;em- 

We will expand the A/ = mass differences assuming that 

1 1 

N Np A ^ ' 

The predictions we obtain will be viewed as predictions for AS ~ 1 mass differences, which 
are 0(e) at leading order in this expansion. 

For the A/ ^ mass differences, we assume that 

nid — m., 

«EM < \ • (2.4) 
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To see that this is reasonable, we note that (to first order in the expansion performed below) 



OEM 



{rud - mu)/A 



30%, 



m^± — mlo 

2 2 

Ms+ - 2M20 + Ms- 



My 



My 



20%. 



(2.5) 
(2.6) 



We therefore expand AI 7^ mass differences to order (m^ — mu)e and oem- Except for the 
fact that we treat as a large parameter, this expansion is identical to the second-order 
chiral expansion that is usually adopted in chiral perturbation theory |jl2,13 . 



D. Effective lagrangian 



The effective lagrangian we use is described in ref. P], and we will not review it in 
detail here. We keep track of the large- group theory by writing the baryon fields \B) as 
elements in a spin-flavor Fock space. The operators that couple to these fields are written as 
r-body operators in spin-flavor space with large- (and Np) suppression factors l/A^'""*"*"^, 
where t is the number of flavor traces used to write the operator [^,0 . The fields for the light 
pseudoscalar mesons 11 are collected in the usual combination ^ = e^^^'^ , where / ~ 114 MeV 
is the kaon decay constant. (Note that the r]' mass is of order Np/N, and is therefore not 
included as a light field.) Explicit flavor breaking appears through the quark masses and 
electromagnetic charge matrix 



rria 



rridl 



N, 



\ 



(2.7) 



E. Power counting for chiral loops 

In this subsection, we discuss the power counting of loop graphs in the expansion de- 
scribed above. This power counting is more complicated than for ordinary chiral perturba- 
tion theory, and it is essential in order to ensure that we have not omitted any important 
contributions at the order we are working. (For example, we will see that there are 2-loop 
graphs that potentially contribute to the baryon mass differences at order rriglnms, and we 
resort to an explicit calculation to see that it does not occur.) We begin by reviewing the 
power counting for "pure" baryon chiral perturbation theory (including only the lowest-lying 
baryon octet). We then discuss the new features that are present in the 1/N expansion. 

In "pure" baryon chiral perturbation theory, a generic term in the effective lagrangian 
looks hke [0] 

where B is the baryon field, and 11 is the meson field. (For brevity, we count oemQ^ as a 
power of iTLq in this subsection.) A loop diagram will contribute to the baryon mass 
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~ (iJvj)' (t)' mM/m,) (2.9) 

where 

C = L + l + {\Du + Sn- Vn) + 1{Db + Sb - Vb) + ^Sb- (2.11) 

Here, L is the number of loops, Dn (Db) is the total number of derivatives in meson (baryon) 
vertices, Su (Sb) is the total number of powers of summed over the meson (baryon) 
vertices, Vn (Vb) is the total number of pion (baryon) vertices, and F{AM/mu) is a function 
of the baryon mass differences and pion masses that appear in the loop diagram. Because 
AM/mn ~ rny^, F can be expanded in a power series in 771^^. The chiral suppression 
factor C is written in this form because Db + Sb — Vb and \Dyi + — Vn measure the 
total number number of "extra" insertions of derivatives and/or powers of iriq in baryon 
and meson vertices, respectively. This is because each baryon vertex contains at least one 
derivative or one power of rrig, and each each pion vertex contains at least two derivatives 
or one power of irtq. 

As an application of this formula, note that at order In we must include 1-loop 
graphs with a single insertion of a 2-derivative baryon vertex (L = \, Sb = 0, and Db + Sb — 
Vb = 1). Since the 2-derivative meson vertices are not measured, we cannot compute the 

In rrig logarithmic corrections to the baryon masses in baryon chiral perturbation theory 

g. 

We now turn to loop graphs in large- baryon chiral perturbation theory. The chiral 
suppressions are as given above, and so we will concentrate on the A^-dependence of diagrams. 
An arbitrary diagram naively gives a contribution to the baryon mass 

^M<A (y^)' (^)^iV^-F(AM/mn). (2.12) 

Because a meson vertex can change the baryon flavor quantum numbers only by order 1 
(rather than A^), AM can have contributions that are at most 1/A^ or rrig. In our expansion 
F(AM/mn) can therefore be expanded in a power series in 1/A^ and m^^. 

Since / ~ N^^"^, eq. (|2.12|) apparently violates large- A^ consistency if Vb — L > 1. 
However, there are cancellations among graphs that make the contributions to the physical 
mass differences at most order A^. Actually, it has never been proven that the required 
cancellations occur to all orders, but many explicit calculations have been done that confirm 
this assertion The 2-loop calculation done in the appendix of this paper provides an 

additional highly nontrivial check.0 



^In fact, the first few terms in the expansion for a given matrix element can be derived by 
demanding that these cancellations occur plH] . 



2-loop corrections to the axial currents are considered in ref. |17|. 
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In our expansion, the argument of F is 



AM 
mn 



-(- 



1/2 



(2.13) 



and we can expand F in a power series in AM/mu for loops of K^s and 77's. (Loops with 
vr intermediate states are neghgible at the order we are working.) The terms with extra 
powers of 1/A^ arising from the expansion of F require fewer cancellations for consistency, 
and we obtain 



SM<A 



167r2/2 



A 



c 



jYmin{Vs-L,l} , j\^m'm{VB-L-l,l} 



A 



rrir, 



1/2 



(2.14) 



The leading contribution to the terms in the square brackets proportional to m^' arise from 
diagrams with n insertions of the 1/A^-suppressed mass term Aq. The largest term in the 
brackets in eq. ( p.l4|) is 



\ rrin 



(Vfl-L-l)/2 



(2.15) 



so that 



'm \C-{Vb-L-1)/2 

5Mb <NK[^ 



(2.16) 



Using Vb — L < L, we can write 

C-\{Vb-L-1)> \{L + 1) + {\Dn + Sn- Vu) + \{Db + Sb- Vb) + \Sb. (2.17) 



Using these results, we can enumerate all of the graphs that we need to consider to expand 
the baryon masses to second order in the expansion described above. From eq. ( p.l7| ) we 
see that we need to consider at most 2-loop graphs. Starting from the minimum number of 
insertions of Aq, it is easy to see that we need to consider the 1-loop graphs of fig. 2, the 
1-loop graphs fig. 3 with zero or one insertion of Aq, and the 2-loop graphs of fig. 4 with 
one insertion of Aq. In the appendix it is shown that the leading contribution of the 2-loop 
graphs cancel, so there is no 2-loop contribution at this order. 



III. EXPANSION OF BARYON MASSES 



Using the formalism discussed above, we now turn to the expansion of the baryon masses. 



A. Leading order 



We first expand the baryon mass differences to order e for the A/ = mass differences 
and to order oem and (m^ — my) for the A/ 7^ mass differences. At this order, there are 
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no chiral loop corrections, and the expansion is determined by the following terms in the 
effective lagrangian: 



6C = -iB\[Ao + A, 



\B), 



(3.1) 



where 



Ai = ai{m} + 



oemA r 
An 



bi{Q}{Q} + b2{Qcr'}{Qa^} 



Here, 



m 



^ I {em,^ + h.c.) , {eQ,e + h-c, 



(3.2) 
(3.3) 

(3.4) 



Note that the coefficients of the operators {Q}{Q} and {Qcr-'}{Qo--'} appear to violate the 
large- counting rule discussed above. The reason is that at the quark level these operators 
arise from electromagnetic diagrams such as the one in fig. 5, which are not suppressed by 
a factor of from gluon vertices. However, in order to obtain a good large- iV limit, 
we must demand that oem ^ '^/N so that the electromagnetic Coulomb energy does not 
overwhelm the strong binding energy for large A^. This also explains why the operator {Q^} 
is not considered at this order, since the coefficient of this operator is the same order as 
the electromagnetic operators considered above.0 At this order we obtain the AI = mass 
relations 



-(S-Ar) + |(S-A) =0, (8%) 

(H* _ s*) - (S* - A) = 0, (3%) 

{Q - E*) - (S* - S*) = 0, (7%) 

(S - A^) - (A - A^) = 0, (40%) 

(S* - A) - (A - AT) = 0, (15%) 



(3.5) 
(3.6) 
(3.7) 
(3.8) 
(3.9) 



and the AI ^ mass relations 



^In fact, this is a good example of how the expansion can differ qualitatively from the 
nonrelativistic quark model. In the quark model, the operators {Q}{Q} and {Q^} are expected to 
be of the same order, while the spin-symmetry violating operator {Qcr^}{Qcr^} would be suppressed 
by I/uiq, where uiq is the constituent quark mass. In the constituent quark model, we obtain the 
additional relation S*~ — S*"*" = S~ — E"*" with accuracy 40 it 15%. Unfortunately, this is neither 
sufficiently successful nor sufficiently unsuccessful to decide whether the quark model works better 
than the 1/A^ expansion discussed here. 
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( 


s--s°)-(s--s+ 


) + (n 




= 0, 


(7 ±11%) 


(3.10) 




A++ - 3A+ 4 


-3A°- 


- A~ 


= 0, 




(3.11) 






) + (n 




= 0, 


(2 ± 22%) 


(3.12) 




2S*° + S*") - (S+ - 


2S° + 


S-) 


= 0, 


(40 ± 100%) 


(3.13) 




(A° - A+ 


)-(n 




= 0, 




(3.14) 




^)-(S+-2S° + S-) 


- 2(n 


-p) 


= 0, 


(3 ± 7%) 


(3.15) 




- A++) + (S*- - 


- 4(n 


-p) 


= 0, 


(40 ± 20%) 


(3.16) 



In addition, we can extract the quark mass ratio 

R=— ^ " =6-— ^ -— - = 110 ±30. 3.17 

md-niu (S+-E-)-4(E*+-S*-) ^ ' 



Here, the name of a baryon denotes its mass. Among the AJ = mass relations, eq. ( p.5|) 
(the Gell-Mann-Okubo relation) and eqs. ( |3.6| ) and ( p.7|) (the decuplet equal spacing rules) 
are valid at order m<j independently of the 1/A^ expansion, while eqs. ( p. 81) and (|3.9|) are 
consequences of the 1/A^ expansion. (In this limit, the mass of a baryon with strangeness 
—S is proportional to S.) Among the A/ 7^ mass relations, eq. (|3.10| ) (the Coleman- 
Glashow relation) is valid up to corrections of order {ma — 171^)171^ independently of the 1/A^ 
expansion [Q. Also, the only corrections to eq. (|3.11|) are second order in isospin breaking 



independently of the 1/A^ expansion. The remaining A/ 7^ relations eqs. ( ^.12| )-( [3.16| ) as 
well as eq. ( |3.17|) are consequences of the 1/A^ expansion. 

The relations are written as linear combinations of mass differences that go to zero as 
iris (for the A/ = relations) or {rrid—mu) 0, Oem (for the A/ 7^ relations). The 
accuracy quoted for the relations (where data is available) is defined as the deviation from 
zero divided by the average of the absolute value of the mass differences that appear in the 
equation. Defined in this way, all of these relations are naively expected to have corrections 
of order e ~ 30%, except for eq. ( ^.llj ), which is essentially exact in our expansion. Note that 



for the A/ 7^ mass differences the relations that hold as a consequence of flavor symmetry 
alone work better than those that depend on the 1/A^ expansion, but there is apparently no 
such pattern for the A/ 7^ differences. 

The value for R obtained is far from the value i? ~ 25 obtained from an analysis of the 



light pseudoscalar masses |T^. The quoted error only takes into account the experimental 
uncertainty of the masses and does not include the theoretical uncertainty from higher-order 
corrections. In view of the large errors in the relations above, we do not take this value very 
seriously. 

A more objective measure of how well these relations work is obtained by fitting the 
mass differences to the parameters given above. The average deviation in a the best fit is 
30 MeV for the A/ 7^ mass differences and 0.35 MeV for the A/ 7^ differences. (If we 
omit the model-dependent value for A" — A++, we get an average deviation of 0.20 MeV for 
the A/ 7^ mass differences.) The fit also gives R ~ 90. 



9 



B. Chiral loops 



The largest corrections to the leading order results in the expansion we are performing 
come from the loop diagrams in fig. 3 and are of order Nm^J"^ for the A/ = mass differences 
and order Nmy^ {md—rriu) for the A/ 7^ mass differences. These diagrams can be evaluated 
from 



5Mb = {B\{TAo']{TAa''}\B] 



pjpk 



6/27 (27r)4 {p2-m\)po' 



(3.18) 



where A = K,r], etc. labels the pseudoscalar mass eigenstates and T4 are the corresponding 
generators normalized so that tr(T4TB) = 6ab- The effects of tt^-t] mixing are incorporated 
by using the generator 



2R 



T^ + 0{l/B?) 



(3.19) 



for the 7] mass eigenstate, where R is the quark mass ratio defined in eq. ( p. 17 ). The result 
can be written 



5Mb = {B\JaOa\B), Ja 



9^ 



where 



Ok 

Or, 



-{N + Nf- 2N,){S} + ^iSa^Ha^} + {S}{S} - ^{Sa^}{Sa^}, 



-{Sa^}{a^} - 



-{Sa^}{Sa^} 



1 

6R 



■{Tsa^}{a^} + {T,a^}{Sa^} 



Ok± 
Here, 



O 



/o 



■Y^r,}-^{n}{S}+^-{na^}{Sa^}. 



S = 



On, 



V 



livj 



\ 



Ons) 



(3.20) 

(3.21) 

(3.22) 
(3.23) 



(3.24) 



We do not need expressions for the pion loops, since they are suppressed by ~ {rriu^/rnsY^'^ 
for A/ = quantities, and by ~ aEMA/m^ ^ for A/ 7^ quantities (since the contribution 
of pion loops is proportional to 771^^+ — 171^^0, which is purely electromagnetic at first order in 
chiral perturbation theory). 

Note that eqs. ( |3.21| ) through ( p.23| ) are valid for arbitrary A''^^^^^ and Np, and that they 
have a sensible limit as A^, Np 00 independently of the extrapolation of the quantities 
Nu,d,s- According to the rules of our expansion, we must keep the full dependence on Nu^d,s 
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and Np since there are limits where each of the above terms is important. The physical 
results are obtained by simply setting Np = 3, N^^a^s = 1-0 

Substituting these expressions into the lowest order relations, we obtain the modified 
A/ = relations 

(S - S) - (S - iV) + |(S - A) = IJk - Jr, 

= -2 ± 1 MeV, (3.25) 
{E* - E*) - (S* - A) = IJk - J, 

= -2 ± 1 MeV, (3.26) 
{n - E*) - {E* - S*) = IJk - J, 

= -2 ± 1 MeV, (3.27) 
(S-iV)-(A-iV) = -|(J^ + J,) 

= -400 ± 160 MeV (3.28) 
(S*-A)-(A-Ar) = |(J^ + J,) 

= 200 ± 80 MeV (3.29) 

and the modified A/ 7^ relations 

(S*- - S*+) - (S+ - 2S° + S") - 2{n - p) = -UJro - Jk+) + 7^^, 



'24\ 

;-4.0±1.6) + (4.8±1.9) (— (0.1 ±0.6 MeV) (3.30) 



RJ 

(A° - A++) + (S*- - S*+) - A{n -p) = -^{Jko - Jk+) + ^J, 

'24\ 



= (-4.0 ±1.6) + (4.8 ±1.9) (^—j (1.6 ±0.85 MeV) (3.31) 

where the errors in the theoretical prediction are obtained by assigning a 20% uncertainty 
to the coupling g = 0.8 extracted from a lowest-order fit to the AS* = 1 semileptonic baryon 
decays, and the experimental values are shown in parenthesis. (The remaining A/ 7^ 
relations do not receive corrections at this order.) In the A/ = relations, the 0(m^/^) 
corrections to the relations that follow from SU{Np) fiavor symmetry alone are small due 
to an "accidental" cancellation. (This was noted in ref. [0.) However, even taking into 
account the theoretical errors, the corrections to the relations that follow from the large- iV 
expansion are too large and have the wrong sign. The same conclusion does not appear to 
hold for the corrections to the A/ 7^ relations, although the situation is obscured by the 
large uncertainties involved. This situation is analogous to what happens for the magnetic 



^Note that for Np > 3 there are "extra" pseudoscalar mesons transforming in the fundamental 
representation of SU (Ng) that also contribute in the chiral loops. The masses of these states are 
determined in terms of the "physical" states by SU{Np) flavor symmetry, and so these contributions 
are calculable. However, these contributions go to zero in the physical limit Np = 3, and so there 
is no reason to compute them explicitly. 
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moments, where the O(my^) corrections to the large- relations are too large, although of 
the right sign [Q. 

One can view this situation in several ways. The most conservative view is that these 
results show that the expansion we are performing does not work well. We will instead 
adopt the attitude that the apparently large corrections at this order of the expansion may 
be misleading, and go on to higher orders of the expansion. At higher orders, these large 
loop contributions can be cancelled by counterterms, and we will see that the predictions 
work very well. 

To see why our viewpoint may be reasonable, we note that there are several calculable 
higher-order effects that we have omitted in our calculation, all of which substantially sup- 
press the correction. The first of these are higher-order corrections to the meson coupling g. 
These are known to be large and apparently reduce g . Second, we have checked that 

using the exact kinematics for the particles in the loops decreases the loop corrections by ap- 
proximately a factor of | for the large- A^ relations. Third, we expect that the meson-baryon 
coupling decreases at large momentum transfer, reducing the effects of K and rj loops. All 
of these effects are suppressed by additional powers of rris, and are therefore higher order in 
the expansion we are performing. Although none of these corrections changes the sign of the 
loop corrections, it is possible that these suppressions are large enough that the expansion 
we are performing makes sense at higher orders. We certainly cannot prove conclusively that 
this point of view is correct, but we can obtain evidence for it by going to higher orders. 

There are also nonanalytic corrections of order m^lnm^. These are formally larger than 
the analytic corrections due to the counterterms by lnA^/M|^ = 1.4 for A = 1 GeV. While 
these effects are definitely enhanced for sufficiently small, we believe that for the physical 
value of rris this enhancement is not numerically large enough to give reliable predictions 
without including the analytic counterterms. 



C. Second order 

In order to expand the A/ = mass differences to order and the A/ ^ mass 
differences to order oem and (m^i — rnu)e, we must include the additional terms 

A2 = ^{ma^}W^} + ^{m'} + W + ^{ma^Hma^. (3.32) 

In addition to the analytic contributions, there are nonanalytic contributions from loop 
graphs of order mglnniq and {mqlnmg)/N. There are m^lnrng contributions from graphs 
such as fig. 2 with the baryon-meson vertices coming from higher-order operators such as 

{{A^A^. ^mm, ^{AVn{AVn. (3.33) 

There are also In rriq contributions from graphs such as fig. 3 with one of the baryon- 
meson vertices coming from the higher-order operators 

i{(VoA^V^}, i{(V,A°)a^}. (3.34) 
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The coupling constants associated with these operators are not well-measured, and we will 
treat them as free parameters. These graphs give contributions to the baryons masses of the 
form 



5Mb = {B\C,0,aKa\B), Ka 



1 



m\ In 



2 ' 



(3.35) 



where the Cj are unknown constants and the independent operators that contribute at this 
order are 



O 



lA 



{Ta]{Ta]. O2A = {TAa'}{TAa'}. 



(3.36) 



There are also {'mq\n.mq)/N corrections arising from graphs such as fig. 3 with a single 
insertion of the flavor- independent baryon mass operator Aq = ^{a^}{a^} /N . These con- 
tributions are proportional to 



[{r^a^},[{r^a^},Ao]]m^ln^ 



m. 



oc 



9' 



167r2/^ N 



{T^a^}{rV} - {TATAa^Ha^}] m^ln^ 



rUq In rUq. 



(3.37) 
(3.38) 



The normalization of these contributions is calculable, but will not be needed. 

When we expand the operators above in terms of the S'?7(iVp') -violating spurious S and 
r3, we see that the operators that appear in Oi and O2 are linear combinations of the 
operators that appear in the tree terms. However, the loop contributions are important 
because they eliminate a prediction for the quark mass ratio R that would otherwise exist 
at this order. 

Eliminating the unknown constants leads to the relations 



(S* - S*) 



(f] - S*) - 2(2* - S*) - 


t-(s*- 


-A) 


= 0, 


(1%) 


(3.39) 


(S* - A) - (S - S) + (S - iV) - 




-A) 


= 0, 


(6%) 


(3.40) 


(H* - S*) 




-s) 


= 0, 


(17%) 


(3.41) 


(S- - H°) - (S- - S+ 


) + (n 


-P) 


= 0, 


(7 ±11%) 


(3.42) 


A++ - 3A+ + 


3A0- 


- A- 


= 0, 




(3.43) 


^^*— '^*o^ (S*^ S*^ 


) + {n 


-p) 


= 0, 


(2 ± 22%) 


(3.44) 


(S*+ - 2S*° + S*") - (S+ - 


2S° + 


S-) 


= 0, 


(40 ± 100%) 


(3.45) 


(A° - A+ 


)-{n 


-P) 


= 0, 




(3.46) 


(AO - A++) + (E+ - 2TP + S') 


- 2{n 


-P) 


= 0, 


(70 ± 30%) 


(3.47) 



Of the AJ = relations, eq. (|3.39|) is an improved version of the equal spacing rule that 
holds to order ml independently of the 1/N expansion |TB| while eqs. ( ^.401 ) and ( |3.41| ) are 
non-trivial predictions of the 1/A^ expansion. These relations were derived in this expansion 
in ref. [0]; the same relations are derived in a different expansion in ref. [Q. Of the A/ 7^ 
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relations, eqs. (p.42|) through (|3.46|) are identical to the lowest-order relations eqs. (p.lO|) 
through ( p.l4| ). It is worth noting that the Coleman-Glashow relation eq. ( 3.42 ) receives 
calculable analytic corrections of order {ma — mu)fns in chiral perturbation theory |]T^; the 
results above show that these corrections are suppressed by 

A fit to the measured mass differences gives an average deviation of 7 MeV for the 
A/ = mass differences (compared to 29 MeV at lowest order). For the A/ ^ mass 
differences, the average deviation is 0.35 MeV, or 0.20 MeV if we exclude the A°-A++ mass 
difference. Since the average of the A/ = mass differences is approximately 150 MeV and 
the average of the AJ ^ mass differences is approximately 3.5 MeV, this is consistent with 
the expected accuracy of l/N"^ ~ 10%. Note that the result for the AJ 7^ mass difference 
is the same as at first order, but this simply refiects the fact that the first-order relations 
work better than expected. 



D. Higher orders 

At next order in large- chiral perturbation theory, we must include the operators 
A3 = ^{m}{a^}{a^} + -^{m} {ma^ } {a^ } , 

which give contributions to the mass differences of order rriq/N'^, m'^/N, and m^. Including 
these terms in addition to the second-order effects discussed above, there is one surviving 
relation: 

(S*+-2S*° + S*") = (S+-2S° + S"). (40 ±100%) (3.49) 

This relation gets corrections from the term 

oemA 1 



Att NA 



{Q}{Qa'}{ma^} (3.50) 



at order a-Eui^s, and also possibly from loop effects at order m^^^, m^lnmg, etc. (Checking 
this requires us to consider 3-loop graphs.) In view of the large experimental uncertainty in 
this relation, we will not pursue the expansion beyond this point. 



IV. CONCLUSIONS 



We have analyzed the baryon mass differences in large-A^ baryon chiral perturbation 
theory with particular emphasis on the chiral loop corrections. One result of our work is 
that the nonanalytic m^/^ corrections to the AJ = large-A^ mass relations appear to be 
too large and have the wrong sign to explain the corrections to the lowest-order relations. 
Nonetheless, accurate results are obtained at higher order in this expansion. At second 
order, there are eight non-trivial predictions of the 1/A^ expansion (which do not follow 
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from flavor symmetry alone). These relations work to better than 10% accuracy, consistent 
with the assumption that the errors are ~ ~ ml/ . 

Ref. |T0| also analyzes the baryon masses in a combined expansion in 1/A^ and flavor 
breaking. The main differences from the present work is that ref. |l^ performs a different 
expansion in which only baryons with strangeness of order 1 in the large- iV limit are con- 



sidered, and nonanalytic chiral loop corrections are not included, 
expansion of ref. [10] are considered in ref. [M. 



Loop corrections in the 



APPENDIX: BARYON MASSES TO TWO LOOPS 

In this Appendix, we compute the 2-loop contributions to the baryon masses shown in 
fig. 4. (The other two-loop graphs can be seen to be negligible using the power-counting 
arguments of section IIIB.) We will find that these contributions are negligible, but we do 
not know any simpler way to see this than by computing them explicitly. 



General Formula 

We begin by deriving the general formula for the 2-loop contribution to the masses in 
the presence of mixing. The baryon self-energy can be viewed as an operator T{E) in the 
spin-flavor Fock space, where E is the energy of the baryon. If we denote the physical baryon 
fields by |n), the mass eigenvalues are determined by 

r(E„)|n) = 0. (4.1) 

This is a nonlinear eigenvalue equation to be solved simultaneously for E^ and \n) order 
by order in the loop expansion. The first step is to expand all quantities in the number of 
loops: 

T{Er,) = E„ - A - Li(K) - L2(K) + ■ ■ ■ , (4.2) 
E^ = + + + ■ ■ ■ , (4.3) 
\n) = |n(°)) + |nW) + |n(2)) + ---. (4.4) 

The energy dependence of the loop perturbations must also be expanded: 

L,(E„) = L,(E(°)) + mE^M'^ + .... (4.5) 

Writing out eq. ( [l.l] ) and equating terms at the same order in the loop expansion gives at 
tree level 

A|n(°)) = Ef |n(°)), (4.6) 

at one loop| 



We are assuming that all degeneracies are lifted by the tree-level terms, so that we do not need 
degenerate perturbation theory. 
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£)(') = (,!<»'|L,(Bf )|«<°)), (4.7) 
|„,.,,^^„„,o,,(l«)^, ,,,, 



and at two loops 



E^^ = + E«(n(°)|L;(Ef )|n 



|(mW|Li(EW)|nW)| 



2 



^ ^ p(o) p(0) • ^^-^^ 

The "energy denominator" terms in the two-loop formula are irrelevant in the limit where 
we can ignore mixing, 

(m(°)|Li(£;(°))|n(°)) = for m 7^ n. (4.10) 

For the problem at hand, mixing always violates isospin, so that the energy denominator 
terms in eq. ( [4.9| ) are second order in isospin violation. This is true for arbitrary and Np, as 
can be seen by considering the generalized isospin and strangeness discussed in section JIB as 
good quantum numbers. Since the baryons of given spin J form an irreducible representation 
of SU{Nf), and because any state in such a representation is uniquely specified by these 
quantum numbers, we see that mixing always violates generalized isospin. We can then 
write 

= (n(°)|L,(Ef) + i [L,(E(°)),L;(Ef )] Jn(°)). (4.11) 

This formula gives the precise meaning of the graphs in fig. 4. The anticommutator term 
can be thought of as arising from wavefunction renormalization. 



2-loop Contribution 

We begin by considering the "true" 2-loop graphs; we will consider the 1-loop graphs with 
counterterm insertions in the following subsection. By the power-counting considerations of 
subsection IIIB, we need only consider graphs of the form in fig. 4 with a single insertion 
of Ao.[| These integrals are regulated using e.g. dimensional regularization so that all the 
manipulations below are well-defined. This contribution is a sum of terms with momentum 
integrals of the form 

fi^. (Po.*) oc (4.12) 

by 3-dimensional rotational invariance. We then obtain 



^We have also evaluated the 2-loop graphs with no mass insertion and verified that the depen- 
dence is consistent. 
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1 



{AXBBA + ABBXA) 



PoiPo + 9o) 
+^r^ ^ABXBA 

+ „ , ^ -AXBAB 

pmipo + qo) 

+ -ABXAB 

Poqoipo + qor 

+ ^-^ -ABAXB 

Poqo{Po + qo) 

- ' + ^ I {AX ABB + BBAXA) 



(4.13) 



V2po?o pm, 

where we have used the abbreviations 

A = {TAa'}, B = {TBa''}, (4.14) 

and e.g. 

ABXBA = ABAqBA - i {AqABBA + ABBAAq) (4.15) 

= l[AB, Ao]BA + lAB[Ao, BA] , (4. 16) 

where Aq = I2{a^ } {a^ } / N . 

Eq. ( [4.13| ) is naively a 5-body operator (if we take into account the commutator structure 
defined by X insertions), but it must be a 4-body operator by large- consistency. By 
explicit calculation, we find that the 5-body part of eq. ( [4.13 ) indeed vanishes, and we 
obtain 

5Mb y / ^ ^ I H H 



36/4 J (27r)4 p2^mlJ (27r)4 g2 _ 
^ ;AB[A,[B,A]] 



Mql 

4go + 3po 
2po?o(po + qo 
2pl + 3p2go + 2pogo + 2go 



A[A,B][B,A] 



Poql{po + qo) 

-|- lower-body operators. (4-17) 

Evaluating the 4-body operators that appear, we find that they are all < iV^, and we have 
checked that all lower-body operators that appear are also < A^^ even for Np ~ A^. This 
shows that 

5Mb < Ar°mjlnm„ (4.18) 
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which is neghgible in our expansion. 

Note that there is a 4-body operator that could have contributed, namely 

AA[B, [B, A]] ~ j^AA {[B, {a^}][B, {a^}] + {a^}[B, [B, {a^}]] + ■ ■ •) 

~ [N'^ + NpN + ■■■]. (4.19) 

Therefore, the important point of this calculation is that this operator does not appear in 
the evaluation of the 2-loop graphs. 

One loop counterterm graphs 

As long as we are evaluating 1-loop graphs, we can use dimensional regularization and 
minimal subtraction, which allow us to simply drop the divergent parts. When we consider 
2-loop graphs, we must be more careful about the counterterm structure, since we must also 
evaluate the 1-loop graphs with a single insertion of a 1-loop counterterm. 

As already stated above, diagrams with L loops generally scale as A^^ for large A^. This 
means that there are divergences in IPI graphs that cannot be cancelled by counterterms 
with the same A^-dependence as the tree-level lagrangian. Even at 1 loop, there is a divergent 
wavefunction renormalization of order A^, while the tree-level kinetic term is order 1. This 
does not imply that the lagrangian is not closed under renormalization, since we can rescale 
the fields so that the counterterms have the correct form.0 We therefore write the (bare) 
effective lagrangian as 

C = {B\Z^'^ . V - A + ■ a} + ■ ■ ■ + at] Z^'^\B), (4.20) 

where Z and 0^ are divergent counterterms to be chosen order-by-order in the loop ex- 
pansion to render all graphs finite. By allowing Z to have arbitrary A^ dependence, we 
can cancel the divergences with counterterms Oct that have the same form as the operators 
appearing in the tree-level lagrangian. This is the sense in which the bare lagrangian has 
the same form as the tree-level lagrangian. 
At one loop, we find 

Z^^--^{1'a''']{Tao'] + 0{A^) (4.21) 

O. = -jf^r..'}. \MA - ,{A ■ .}]] d.v /|^ -L. (4.22) 

Here, "div" indicates the divergent part. 

The counterterms in Oct are smaller by one power of 1/A^ than required by large- A^ 
consistency. This can be used to see that the contribution of the 1-loop counterterm graphs 



^Again, this has not been proven to all loops. 
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to the baryon masses are ~ iV" (not ~ N). The reason is simply that the physical mass is 
independent of the scale of the fields, and so one can compute the masses with Z = 1. All the 
other counterterms that appear are now suppressed by an additional power of 1/N, which 
is immediately gives the result that these contributions are negligible in our expansion. 
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FIGURE CAPTIONS 



FIG. 1: The Young tableaux for the SU{Nf) flavor representation of the spin- J baryon 
multiplet. 

FIG. 2: A contribution to the baryon mass. The solid lines are baryons and the dashed lines 
are mesons. 

FIG. 3: A contribution to the baryon mass. 

FIG. 4: Two-loop contributions to the baryon mass. 

FIG. 5: Example of a quark graph giving rise to the 2-body operators ^-iid 
{Qcr-'}{Qcr^} in the effective lagrangian. The wavy line represents a photon and the curly 
line represents a gluon. Note that there is a factor of N from the color sum, so this graph 
is order 1 in the large- iV limit. 
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